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ABSTRACT 

This report covers work done on performance analysis of a 5 inch, 
54 caliber gun-launched guided projectile with solid fuel ramjet (SFRJ). 

A computer program (TRAJET) was developed. The program contains 
ramjet and trajectory analysis. The ramjet part considers conical 
shock wave losses, inlet boundary layer losses, normal shock losses, 
subsonic diffuser recovery, expansion into combustor losses, heat losses 
at the combustor and nozzle losses. 

A flat earth trajectory with drag and thrust was considered. The 
various drag coefficients which were considered are: cowl drag coefficient, 
skin drag coefficient, wing (or fin) wave drag coefficient and wing (or 
fin) friction drag coefficient. Base drag is assumed to be zero due 
to the jet from ramjet nozzle. 

It was found that the cowl drag coefficient has a major influence 
on the results. Therefore, a separate program (AERO) was developed to 
calculate this parameter. 

The 5"/54 solid fuel ramjet has a capability to produce fuel 
specific impulse in the order of 400-900 sec. depending mostly on the 
flight altitude. The thrust coefficient varies in the range of 0.3+0.1 
depending on the internal areas. 

A range in the order of 50 miles can be achieved with the ramjet 
operation compared to only 13 miles achieved by the conventional pro- 
jectile. At low-altitude launch, a range of over 18 miles can be reached 
with the ramjet version. Launches at nigh elevation angles can be 
useful in air-defense scenario. The ramjet propelled projectile reaches 


the ranges mentioned above at high Mach numbers (My>1.8). It is, therefore, 


clear that the ramjet concept provides significant improvement and has 


an Anti Ship Missile Defense (ASMD) capability. 
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1. INTRODUCTION 
1.1 Backaround 

This report covers work done on propulsion and flight me- 
chanics of a gun-launched guided projectile. Gun-launched guided pro- 
jectiles have been developed by Martia Marietta for the U.S. Army 
(Copperhead) [1,2] and more recently, also, for the U.S. Navy [3,4] 
(5"/54 Mark 46). The USN round has solid rocket propulsion. 

The addition of a liquid fuel ramjet (LFRJ) to the Navy's 
version, was examined in the past by Brown [5]. This report concentrates 
on the addition of a solid fuel ramjet (SFRJ), instead. 

It is believed that solid fuel ramjet has some potential 
advantages compared with the liquid fuel. Some of these advantages 
are: 

-Simple desion 

-High reliability in operation 

-Low cost 

-Fuel control system not needed 
On the other hand, there are also a few disadvantages to SFRJ compared 
to LFRJ. Some of these are: 

-Difficult to control magnitude of thrust 

-Difficulties in achieving high combustion efficiencies 

In both cases, the addition of propulsion imoroves dramatically, 
the performance of the projectile by multiplication of range and enhanced 
maneuverability. Even more; to operate and produce thrust, the ramjet 
engine depends only on its forward motion at supersonic speeds and does 
not employ any moving parts. This fact, which is especially emphasized 


in SFRJ, leads to some advantages of the ramjet concept over the other 


Propulsion alternatives, at supersonic speed. On the other hand, the 
ramjet engine requires an auxiliary booster to accelerate it to its 
Supersonic operating regime. The boost required causes some system 
difficulties. But, while solving these problems, the ramjet system 
becomes even more attractive for use with gun-launched guided pro- 
jectile, like the U.S. Navy 5"/54. 

A computer program was developed to analyze the performance 
of the SFRJ. The computer program was written for the IBM-370 
computer at the Naval Postgraduate School, Monterey, California. HTPB 
was selected as a fuel, but performance with any other fuel can be tested, 
uSing the same model. A flat earth trajectory with drag and thrust was 
considered. Using solid fuel, a thrust-equal-drag trajectory is more 
difficult to achieve with SFRJ. Therefore, most of the results given 
in this report eliminate this case, and the exact value for drag, at 
each point, was calculated. However, if desired, it is possible to 
change the air mass flow, in order to obtain thrust-equal-drag flight. 
The computer program can calculate this case also. 

1.2) Basie! Concept 

The ramjet engine consists of an air inlet, which serves as 
a diffuser, a combustion chamber, and an exhaust nozzle [6,7]. The 
diffuser admits air to the engine, which is mixed with fuel (solid or 
liquid) at the combustor. After the burning process, which adds heat 
to the flowing air within the system, the gases are transfered to the 
nozzle. The nozzle converts part of the thermal energy into kinetic 
energy to produce thrust. 

The areas inside the ramjet engine are usually divided into 


six stations, as illustrated in Figure 1.1. Station 0 defines the cross- 
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section area of the stream tube captured by the inlet of the projectile. 
Stations ] and 2 identify the diffuser. Station 1 itself is at the 
throat of the diffuser, but sub-station 1C is after the conical wave; 
sub-station 11 and Sl are ahead of the normal shock wave, located at 
the throat, or at the actual place, respectively; sub-stations 12 and 
S2 are as above, but behind the normal shock wave. Stations 3 and 4 
refer to the entrance and to the exit of the fuel grain within the 
combustor, respectively. Note that both Ay and Aa increase with time 
as fuel burns. Station 5 and 6 belong to the nozzle's throat and to 
the exit of the nozzle, respectively. 
1.3 Data 
1.3.1 Dimensions 
In order to be compatible with the Navy's 5 inch, 54 caliber, 
Mark 46 gun mount, as modified for gun launched guided projectiles 
(Figure 1.2), a set of requirements were adapted initially. These were: 
a. External shape of existing 5" guided SAL projectile 
b. Length - 60.97" 
c. Length of combustion chamber - 23" 
d. Total weight - 104.7 1b. 
e. Muzzle velocity - 2500 ft/sec. 
Typical values for the internal areas in the ramjet within the 
projectile are (units - sq. in.): 

A. Ag Ay Ay Az As Ag 
19.3 5.2 220 4.3 8.2 as 1s 
Refer to Figure 1.1 for definition of the areas. A, 1S a reference area. 
For typical flight Mach number of: My = 3, the appropriate Mach numbers 


at the main stations are typically: 


M M 


Ne ll 12 2 3 5 
Dae Diol 0.56 0.3 O71 ] 2 


1.3.2 Combustion Process 
Some of the losses in the total pressure were taken to be 
constant. These are: 
a. Inlet boundary layer losses (T' ys typically = 0598 ye 


b. Subsonic diffuser recovery (17 93 typically = 0.93). 
c. Nozzle losses (m3 typically = 0.96). 
These values are typical and are not expected to vary too much. 
All the other losses in total pressure, which are dominant to 
the projectile performance,were calculated. These are: 
a. Conical shock wave losses (.). 
b. Normal Shock losses (Tr). 
c. Losses due to expansion into the combustion chamber (m,). 
d. Heat losses in the combustion chamber (m)). 
See section A 1.4.2 for definition of various 7 . Combustion efficiency 
was taken constant: Ny = 0.90 


Air heat capacity ratio was also taken constant: = 1.4; however, the 


Ya 
value for the gas heat capacities ratio of the combustion products (y,) 
‘was calculated from thermodynamic data for Hydroxy Terminated Polybutadiene, 
HTPB, burned in air [8]. The stoichiometric chemical reaction of HTPB 
burning with oxygen is as follows: 

4 Cray 939, + 393 0, + 292 CO, + 206 HQ0 oe Gea) 
1.323 Trajectory 
In the trajectory part of the program, the various drag coef- 


ficients were calculated. Those are: 


a. Cowl drag coefficient (Con) » [9-12]. 


b. Skin drag coefficient (C,.), [13-15]. 


DS 


c. Wing (or fin) wave drag coefficient (C 16 - 20}. 


pwn)? | 

d. Wing (or fin) friction drag coefficient (Cove) s [134 49%. 
Base drag [13, 21-23] is assumed to be negligible due to the jet from ramjet 
nozzle. 

The model, which was chosen to calculate the cowl (nose) drag 
coefficient, was based on a theoretical development done by T. H. Gawain 
[9]. The modified program (AERO) is listed in Appendix G. However, 
program AERO as it is, appears to be too long to be used directly in the 
main program (TRAJET). Therefore, best fit curves for calculated results 
from AERO were used in TRAJET. 

Skin drag coefficients were calculated for either laminar or 
turbulent flows. The same routine was also used to calculate wing or fin 
friction drag coefficients. To calculate the wing wave drag coefficients, 

a psuedo 3-dimensional model was developed. 

The program also has an option to calculate a trajectory of a 
projectile without propulsion. In this case, the drag coefficients which 
are calculated are: 

a. Nose drag coefficient (a different model than the above). 
b. Base drag coefficient 
c. Skin drag coefficient (as above). 

1.3.4 Air Defense Scenario 

In the air defense scenario, the program takes into account only 
cases in which the projectile exceeds a Mach number of at least 1.8. This 
value of minimum Mach number (XMy) can easily be changed. 

1.4 Results 
Each section of the program was developed, tested and run sep- 


arately. The ramjet part was first run without the trajectory part using 


constant altitude (Typically - 10,000 ft). The same was done with the 
trajectory part using vacuum case, thrust-equal-drag flight, or constant 
thrust case. The final version was programmed to give the following 
optional printings: 
a. Loop on all possible values of Ad/A. and Ac/A,, and print 
summary tables only. 
b. Print detailed, time dependent, tables for any specific area 
ratio chosen: 
- Results from combustion process (file name: CMB D) 
- Results from trajectory process (file name: TRU D) 
- Various drag coefficients (file name: DRG D) 
c. Detailed print of every step during the calculation, for 
checkup. 
d. Variation of the above: 
Detailed print of cases that were found not to be 
suitable: 
- Reasons only 
- Full detailed parameters 
- Loop on Mach numbers, also (output of subroutine 


CALCM) 


2, YPICAL RESULTS 

Figure 2.1 presents the dependence of the fuel specific impulse (Teo: 
in sec.) of the ramjet on the projectile Mach number at various altitudes. 
It appears that the 5"/54 ramjet has a capability to produce fuel specific 
impulse in the order of 400 - 900 sec., depending mostly on the flight al- 
titude. The dependence of = on the flight Mach number is weak. 

Figures 2.2 and 2.3 present the dependence of the thrust coefficient 
(C,) on the internal area ratios Ag/A,. and Ac/A,.. In figure 2.2, the change 
of Ce with altitude and with Mach number is also presented. The thrust coef- 
ficient (Ce) varies in the range of 0.3+0.1 while Ao/A,. changes from 0.25 to 
0.40 and Ac/A,, changes from 0.42 to 0.26. 

The correlation between the fuel specific impulse (I.,) and the thrust 
coefficient (Ce) is presented in figures 2.4 and 2.5. In both figures, a Mach 
number of My = 3.0 was selected. In figure 2.4, the correlation was checked 
at various altitudes and at various Ad/As area ratios. In figure 2.5, various 
internal area ratios (Ag/A.» A./A.) are presented. 

More detailed results are presented in Appendix H. The dependence of the 
projectile performance on the other internal area ratios (A,/Ag; Ao/Ags A2/A,) 
was also checked. Some typical results are presented in that Appendix. 

An altitude ys range dependence for various elevation angles is presented 
in figure 2.6. A range of over 80 km can be achieved with the ramjet operation, 
compared to only slightly more than 20 km achieved by the conventional pro- 
jectile. A low-altitude launch (in this figure, an elevation angle of 15° 
was selected) is also presented reaching a range of over 30 km with the ramjet 
operation. The high elevation angles are mostly used in air-defense scenario. 


The drag of the projectile when the ramjet is not operating, for example, 


after burnout, was not determined. The computer program does not account 
for the drag increase due to the ramjet not operating. Consequently some 
of the trajectory curves in figure 2.6 are in error. However, for tra- 
jectories at low gun elevation, the ramjet burns all the way to splash. 
These trajectories are accurate. For all trajectories, the curves are 
accurate to the point of ramjet burnout. The trajectories of interest to 
air defense are accurately calculated. Trajectory of thrust-equal-drag 
(vacuum) case is shown for comparison. 

Results for air defense scenario are presented in figure 2.7. Only 
ranges where the projectile Mach number exceeds at least 1.8 where con- 
Sidered. The area ratios where chosen as specified in the figure. The 
two cowl angles, shown in figure 3.1 were 20° and 9.5° respectively. 

The gun elevation angle was varied from 7° to 80°. The change of atmospheric 
conditions with the altitude was taken into account. In table 2.1, some 


typical results for "Surface-to-Surface Mission" are presented. 
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TABLE 2.1 


Solid Fuel Ramjet Propelled 5"/54 Projectile 
Surface-to-Surface Mission 
Ranges (km) vs Gun Elevation Angles 


Elevation 
Angle rie 25° 45° 65° 80° 
By eamee Be 49.9 80.3 15.9 56 
b. Projectile ONZ 75 20R2 16.3 7.6 
Without 
Propulsion 


Note: 1. Area ratios as in Figure 2.7 
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3. Discussions 

Looking back at figure 2.1, the dependence of Ty on altitude and on 
Mach number should be explained. We shall do that by using the equations 
described in Appendix A. 
From equations: (1.4.9a), (1.4.12), (1.4.16) 
together with equations: (1.1.1), (1.1.2), (1.1.3) 


one obtains: 


Ie, = Ca A Ko (3.1) 
where: X = eis (322) 
and: Ce = Xo - ko (3.3) 
where: et oe (3.4) 
X, = kgMyl1+ke (Py Mg) ] (3.5) 


The parameters kK, to kp are functions of the various area ratios, the 
temperature of air (To); the heat capacity ratio of air Sey, and the 
perfect gas constant (R.). These parameters are assumed to be constants 
in discussing the influence of the change in altitude and in Mach number 
on the value of Ley: The altitude dependence is mainly due to change in 
atmospheric pressure (Pq). In the conditions chosen for figure 2.1, the 
dependence of Ion on pressure is approximately a at My =3. That means 
that, in the region mentioned, ley pressure dependence is mostly due to 
change of X (equations 3.1 & 3.2). From the same equations, the Mach 
number dependence of = can also be explained. At high altitude, the 


1.41) 


change of X (My is very close to the Ce dependence on Mach number 


(m)7' 9) and therefore 5p is almost constant while changing M 


0° 
other hand, at sea level, the change in rer Pag) is smaller than that of 


On the 


-2.14 . 
Ce(My ) and therefore re changes with My as shown in figure 2.1. 
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Figures 2.2 - 2.5 present similar dependences of the ramjet performance, 
and can well be understood using the same equations. 

Testing these results, together with those presented in Appendix H, 
the design of the ramjet internal area ratios can be completed. The 


results are listed in Table 3.1. 


20 


Table 3.1: Ramjet Design 


Dimensions: 


External diameter = 5" 


Total length 60.97" 


Total weight 104.7 1b (47.5 kg) 


Area Ratios: 
A/a. Ay/A Ay/Ag A3/A., 
0.25 0.47 0.887 0.426 0.26 
Reference Area: 
A = 19.3 sq. in. (124.5 om*) 
Combustor 
Solid fuel: Hydroxy Terminated Polybutadiene (HTPB). 
Fuel weight: 3 kg 
Fuel density: 971.56 kg/m> 


): 400 - 900 sec. 


Fuel specific impulse (I, 


Booster 
Booster weight: 2 kg 
“Booster density: 1650 kg/m° 


Booster specific impulse (I..): 240 sec 


sp 
Performance 

Muzzle velocity: 762 m/sec 
Velocity after booster: 863 m/sec 


Thrust Coefficient (Ce): 0.2 - 0.4 
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In figures 3.1] - 3.4, the designed ramjet concept is presented. The 
guidance and the control sections as well as the warhead were not redesigned. 
The location of the tailfin is described in figures 3.3 - 3.4 [White,24]. 
The configuration of the solid fuel ramjet presented here, is in accordance 
with the design of the liquid fuel ramjet done previously by Brown [5]. 

The new design is also in agreement with the Navy's requirement to be 
compatible with its 5"/54 Mark 46 gun mount, as modified for gun launched 
guided projectiles. 

Figure 2.6 presents a full-range comparison of the SFRJ 5"/54 pro- 
jectile performance with that of the conventional projectile. The 
improvement in performance of the ramjet-propelled, guided projectile is 
very significant. The extended range of the ramjet concept can be used in 
a "Surface-to-Surface Mission" (Table 2.1). The main improvement of the 
ramjet concept might be in an "Air-Defense Mission" (figure 2.7). A 
minimum Mach number of My = 1.8 was assumed for both the ramjet concept 
and the conventional projectile. At lower altitudes, the projectile is 
limited by burn-out of the fuel. At higher altitudes, the performance 
boundary is due to ramjet unstart. The conventional projectile is limited 
by the Mach number decay to values less than 1.8. It is self-evident 
that the ramjet concept provides significant improvement, and, therefore, 


has significant ASMD Capability. 
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MAXIMUM OUTSIDE DIAMETER 
<1 ——— 


ALL DIMENSIONS IN INCHES 


ALLOWABLE AREA 
FOR EXHAUST 
NOZZLE 


SECTION A-A 


Figure 3.4 Aft Body Fin Design: Section A-A [24] 
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4. Conclusion 

Based on the computer model discussed above, the ramjet-propelled, 
guided projectile provides significant improvement over the conventional 
projectile. The fuel specific impulse of the SFRJ is in the order of 600 
+200 sec, depending mainly on the altitude of the projectile. The 
appropriate value of the rocket is only 300 sec. The thrust coefficient 
varies from 0.2 to 0.4, depending on the atmospheric conditions (altitude) 
and on the geometry of the projectile (internal areas). Therefore, the 
ramjet-propelled, guided projectile reaches a range of about 80 km compared 
to range of slightly more than 20 km in the conventional projectile. The 
improvement of the ramjet concept might be in both "Surface-to-Surface 
Mission" and in "Air-Defense Mission". It provides an ASMD weapon which 
is complementary to guided missiles. 

More work is required to design an optical system incorporated into 
the inlet. For satisfactory ramjet performance at the flight Mach numbers, 


the lens must be a conical shape. 
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Appendix A: SOLID FUEL RAMJET: EQUATIONS 
Al. Combustion 
Al.1 Computation of fuel - regression rate, weight rate of 
burning fuel, and fuel - air ratio 
Define G as the weight flow rate of air per unit through the 


entrance port to the combustor; see figure 1.1. Hence G is agiven by: 


We 
ole Gana) 


and has dimensions of Ib/sec.in’. The simple form of the regression 
rate of a solid fuel, *, iS given by: 

f = aG" (1.152) 
Where a and n are empirically determined constants. The dimensions for 
fr are in/sec. Knowledge of the temperature dependence of the regression 
rate will allow the use of a more accurate model instead of equation 
(Tee 
The weight of fuel burned per unit of time is as follows: 

We = o¢ trD,L, Chelle) 
here p, is the density of the fuel in 1b/in®, Dj is the inside dia- 
meter of the fuel grain. Note that De increases as the fuel burns. 

The lenath of fuel grain is Ly inches. For the case of HTPB, the value 
of 0, 035i1 1b/in? was taken. In an actual solid fuel ramjet, * varies 
alona the grain; * is largest in the region immediately downstream of 
the entrance port of the combustor. However, for the ramjet model 


developed here, the value of * is assumed to be constant along the grain. 


Consequently: 
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The integral form for the change in grain internal diameter is used 
since * may vary with time. The initial grain inside diameter is Do 
inches, and the inside area for fuel grain as a function of time can be 


written as: 
= rt 
le ee | 
A, = 4 A Bag + 2) rat | (1.1.5) 


By definition, the fuel - air ratio is: 
Ne 
f = aT (1.1.26) 
a 
The value for W. is obtained from weight flow through the inlet, and the 
value for We is calculated using equation (1.1.3). The total mass flow 
through the nozzle is given by: 


We = We ~ We = Ww. + f) Gee ys) 


For HTPB burning in air, the stoichiometric value for f is 0.0728; high 
combustion efficiency is difficult when f is less than 0.025. 

Al.2 Computation of Combustion Exit Condition 

Combustor exit conditions are specified by four quantities as 
follows: stagnation temperature, Tra> On5 stagnation pressure, Pras 
psi; ratio of heat capacities, Yes and gas constant, Res in/Op. In 
the computer program, the appropriate mks units are used, i.e. re kg/m, 
m/O,, respectively. To determine the exit conditions, two input quan- 
tities are needed. These are fuel - air ratio and stagnation temperature 
at the combustor inlet, Tro: Note that Tr3 = Tr9 has been assumed. 

From the thermodynamic data for HTPB burning in air, one determines 
me ( th) 


R The symbol Tyq(th) is a theoretical temperature which 


T4 ’ Tee f° 
results from 100% combustion efficiency. Introducing the definition of 


combustion efficiency yields: 


ipo) 
ive) 


er =} Tra(th) - ina, Der (1.2808 


As discussed previously, a constant value of Nt equal to 0.9 has been 
assumed. 

The value for Pra is calculated based on one-dimensional, choked 
nozzle flow. A certain value of Prq is required to force a certain weight 
flow, He » through the nozzle. Assume that Y¢ remains fixed through the 


nozzle. Define a function of Ye as: 


, \Y¢ a 1)/12(y¢ - 1)] 
T = VV¢ —_= | (la2e2) 


Define a characteristic nozzle velocity, 
c*, m/sec: 

SOTA 
8 seege 


a (1.203) 


where g is the acceleration of gravity and has value of 9.807 m/sec’. 


The required value for Pra is given by: 


(1.258) 


The decrease of flight stagnation pressure, Pro» by inlet and combustor 
losses must not be too large. If the inlet and combustor do not provide 


the required Pra» the inlet will unstart and We will decrease. 


Al.3 Computation of Nozzle Exit Conditions 


The relation between the area ratios and the Mach number’ is well 


known by the formula: 


Ci) 


As. Ag are the areas at the throat and at the exit of the nozzle, re- 
spectively. Knowing Yes the exit Mach number (Me) can be calculated 
for any nozzle area ratio (Ac/Ag). This indirect calculation is done 
in subroutine CALCM, using Newton - Raphson's iteration routine. 

The total pressure at the exit of the nozzle (Pre) is defined by: 


Prg = Prq ™ Ulgsee) 


where the total pressure at the exit of the combustor (Pry) was calculated 
previously (e.g. 1.2.4). 
Knowing the total pressure at the exit of the nozzle (Prg) and the 


Mach number at this point (Me), the exit pressure (D¢) can be calculated: 


Yent 9 Ye) ter!) 
a Pr6 1 ame am (1.3.0) 


A1.4 Computation of Thrust and Thrust Coefficient 
Al.4.1 Thrust Coefficient (C_) 


The thrust of the engine is the net rate of change in 


momentum at a steady state condition, and is given by: 


b= Peg + they - Pog - mhU 


Pot + tee - Porg - Muy 


An - Ag) + Dg (Ag-Ag) 


0) o| 4 6 
(1.4.1) 


9 ~ Pol 
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where Ug» U 6? tho» the are the velocities and inass flow at the inlet 
entrance and at the nozzle exit, respectively. 


From the continuity equation, the following relation arrives: 
2 


mU = UTA Gide) 
Substituting for the density (c) from the perfect gas equation of state: 
5 
RT Claes) 
gives: mU = on ueA (aeae ay) 
From the definition of Mach number and speed of sound: 
US = Meat ; ac = yRT (1.4.5) 
Therefore: nU = BoméyRTA 
fU = pMoA Gee 
Substituting (1.4.7) into (1.4.1) gives: 
A 
S 2 6 2 
F = PgAg (1 + yee ) - PoAg a + y aig ) edes8)) 
The thrust coefficient is defined : 
Cees 
fi doh. (1.4.9) 
A = Oey 2 
where: qo = 0 qUy = EYP gM (1.45 10) 


Combining equation (1.4.9) and 1.4.10) gives: 


CS ae (1.4.11) 


v6 2 
~4 Povo Ay. 


Substituting for the thrust from equation (1.4.8) turns equation (1.4.11) 


into: 


2A-/A 
Ce = oe Fe a + ¥oM Sonja! (iar 
YaMo~ -L Pre/Pe fee A 
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Al.4.2 Pressure Losses 
The pressure ratios in the above formula: 
Pr6/ Pro 


Pr6/P6 
can be substituted by a function of pressure losses across the ramjet: 


Pre _ Pre Pra Pr3 Pr2 Pro (1.4.13) 


Pg Pra Pr3 PT2 Pro Po 


We define the pressure losses as follows: 


p 
ee 7m. = Nozzle losses 
Pr4 
PTA | : 
7 ae Rayleigh flow losses 
P73 
Pr. - 
— = mt, = Combustor expansion losses 
PT? : 
ae _ _ ; ‘ ‘ 
nl Inlet losses = (conical wave loss )*(boundary layer loss) 
TO 
*( normal shock loss )*(subsonic diffuser recovery) = Tet) TycTy" 
Therefore: 


°T6 = lo i Pro = he Pro (1 4 14) 
pee a ee Ns Po i = losses Jpg a. 


Finally, substituting for total pressure ratios (Pr9/ Pg; Pr¢/P¢) gives: 
er at) 


1+ 3M 
Pre/Po _ [tony | i (Yg-1) (1.4, 15) 


Pre/Ps j=losses i, n one 


This relation can be used in equation (1.4.12) which calculates the 
thrust coefficient of the system. 


Al.4.3 Computation of Thrust 


—_—. 


The thrust can easily be calculated from the thrust coefficient, 
using eq. (1.4.9): 


b= C A (1.4.9a) 


f 49 ">, 
The dimensions of F are Newtons (after multiplying eq. (1.4.9a) by 


the acceleration of gravity, g). The fuel specific impulse, Tey» in 


N/kg/sec, is defined by: 
55 = F/We Gla 4-nlGy) 


The specific fuel consumption, SFC, in kg/hour/N is given by: 
SFC = 3600/1. , (liad cili7o) 
Ramjet performance is specified in terms of the performance parameters 


Co, ns be wand’ SEC: 


ihe sp 
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A2. Check for Choked Nozzle 
The total pressure at the throat of the nozzle is given by: 
Prs = Ptq I" (2.1) 


Again, Pra» is the total pressure at the exit of the combustor, and is 


calculated by eq. (1.2.4). am. is the nozzle loss. It follows that: 
2 \¥el (ygtl) (2.2) 
P5 ~ PTS (77) 


The pressure at the throat of the nozzle (pe) should be equal or greater 
than the atmosperic pressure (Dy). - 
P5> Po (2.3) 
When inequality (2.3) is satisfied, the nozzle is choked. 
A3. Heat Losses at the Combustor 
A3.1 Mach Number 

a. Continuity 

At the combustion chamber, fuel is added and the continuity 
equation is: 

03 Uz Ay (1 + F/a) = Oy U, Ay Giisl) 

where As and Ay refer to the entrance and to the exit of the combustor, 
respectively. By assuming that Ay = An» the continuity equation 


(eq. 3.1.1) can be written as follows: 


(ex 


O 
(1 + f/a) ee (33.2) 


“4 3 


Replacing the velocities (U, and Us) by the appropriate Mach numbers 


(eq. 1.4.5), turns eq. (3.1.3) into: 


@) M | i ey. R 
(1 + f/a) des eet . (3.1.3) 
“4 . 34/ '3 Ya “a 
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Where RD Re are the gas constants of air and of the combustion products, 


respectively. Hence: 


a 
("2). (28) ($2) (ver) (_a_ 
Ma P3 T3 VARS Ll + fla (27484) 


b. Momentum 
Applying the conservation law of momentum to the discussed 


problem: 


ee 2 
P3 - 23 U, = Py + Py Un (3s D)) 


From the definition of Mach number (eq. 1.4.5) and the perfect gas 
equation of state (eq. 1.4.3): 


2 


Overt = Ne P; M. (3136) 


Substituting equation (3.1.6) in equation (3.1.5): 


2 
Dob ital 3h Gan 


Pele Ye My 


Again, from the perfect gas equation of state: 


Dies ie (3.1.8) 
Pa. 837 \'3 


Substa tution ef (3.1.8) into (3.1.7): 
2 


(=) 2 (=) eee (3.1.9) 
T 
23 ay) Na vey’ 


Substituting (eq. 3.1.9) turns equation (3.1.4) into: 
2 


“2 
BY CE RE al Rg le 
M) \T3 ee a On| ee T + t/a (3.1.10) 


veMg 


t|K 
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Replacing the temperatures T3 and Ty by the appropriate total temperatures: 


and assuming that the change in total temperature up to the entrance of 


the combustor, is negligible ( , one obtains: 
*3 
% 
rey -\'T0 ‘ree vey ar fie are : 


c. Solution 
The Mach number at the exit of the combustor (My) can be 


solved knowing the conditions at the throat of the nozzle: 


Ae { me (yet1 )/[2(¥~-1)] 
al Scam cca (3.1.12) 


The computation is again indirect, using subroutine CALCM. Knowing 
Mas Equation (3.1.11) is used to compute M.. The solution is received 
by iteration. As first approximation, M. can be solved from the 


following equation: 


i. ( + “at = 8 (3.1. 11a) 
ms 
where: =f 1 + y—My (=) as (1 + #/a) (0) 
B = M 7 Y R 
nl T0 agi 
consequently: rae 
+B- /B -4y 5 
Fk. SAIN Bis 3 A 
Mon = 2, (eu) ie) 


The subscript N shows that the computation was made from the nozzle 
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direction. Now, B is changed to be: 


| Yar! 2 3 
1 + Spe ie 
Besa Sees ieee / (32.1515) 
Ge 2 
fe et | 


This expression is consistent with equation (3.1.11). Substituting 
back into (eq. 3.1.14) gives an improved value for Ma. This procedure 
can be repeated several times, but it was found that even after two 
interations, the value received for Mon is accurate enough, due to the 
smal1 change in B resulting from (eq. 3.1.15). 

A3.2. Total Pressure 


By definition: 


D> : vin! gf (y¢-1) 
Tol il 2) 
DP ) 2 j ( (Sie 
Consequently: 


ete) | 
| Wal pee 
Rw ices uubaira is (3.2.2) 
Prgets /(Ra ey en area 
Ne a 


Substituting for (D4/P,) from equation (3.1.7) results: 


( ‘ 2 | y.-1 ale (y+ ) 
i oan Lert 
Pian Pra go ah 


[? ? | Veet 5 ey (3.2.3) 
Ee 
| 4 


Where the subscript N was defined previously. 
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A4. Computation of Mach Number and of Total Pressure at the Various 
Stations of the Inlet 


A4.1 Initial Conditions 

In the previous sections (Al - A3) the pressure conditions at 
the combustor and at the nozzle region were calculated. Here, the 
pressure conditions at the inlet will be calculated independently. 
Knowing the total pressure conditions at the various stations of the 
inlet will allow the check of whether the inlet can supply the amount 
of air needed by the combustor. As will be seen afterwards, this check 
will also allow to specify the location of the normal shock wave at the 
inlet. 

Assuming that the static pressure (Dp), the static temperature 
(T9) and the flight Mach number (M5) are known from the trajectory part 


of the program, the total pressure and the total temperature can be 


calculated: 
¥,71 y,/(y,-1) 
Prog = Po fle + Moo] : (A200) 
tan 2 


The weight flow through the inlet is given by: 
Wa = Pg Yo Ap (5) (4.1.3) 


Usually, when flight Mach number is equal or greater than the inlet 
design Mach number, the value for Ac/Ag is unity. But for flight Mach 
number less than design Mach number, Ac/Ag becomes less than 1.0. A 
value of 0.9 was selected as a constant value for Ac/Ag- The additive 


drag due to Ac/Ag < 1 was ignored. 


39 


A4.2 Conical Shock Wave Loss 
In this section, the conical shock wave loss will be computed; 
the calculation results include the total pressure, Mach number and 
area behind the conical wave (Price Mies Aro: respectively). 
A4.2.1 Pressure 


The pressure coefficient can be defined as follows: 


Dy --D 
c, = + _, ea (a 
x Po Mo 


For a cone, the pressure coefficient (C,) can approximately be for- 


1.69 
rh 0.096 fo} 
Cc, = | 0.083 a We | (#5) (AeZee) 
0 


where a is the cone half angle. The difference in pressure on surface 


mulated as: 


and behind shock wave is ignored in this model. 
Knowina Cs the pressure ratio (Py ¢/Po) can be calculated (4.2.1) 


ae (zea) 


The same pressure ratio, is also related to the Mach number, normal to 
the conical shock wave (M_). Using (eq. 2.48a) in ref. [16], one can 


get: 
-1) (422.3)) 


—— =] + + 
D9 itl 


Knowing (D1 ¢/Pg) from equation (4.2.1a), M, can be calculated from 


equation (4.2.3) : 


Bie yt 4 (4.2.3a) 
Mos i ice . i aa 
n Po ey. 
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After computing the pressure ratio due to conical shock wave (equation 
4.2.la) and the Mach number normal to the cone (equation 4.2.3a), one 
can use equation (2.54) in reference [16] to compute the total pressure 


ratio at the conical shock wave: 


2 vy AG -1) 
D 2y -1/(y.-1) y.+1, M tar 1 
Te = it = f + = 7 (M 21) | : KE - C5 | (4.2.4) 
“TO fa B rs M +2 
From equation (4.2.4) one obtains: 
Pric = to ee esta) 


A4.2.2 Mach Number Downstream of Conical Shock Wave 
Using figures 4.1 and 4.2, the wave angle (8) can be defined 


as: 


M 
8 = arcsin 4 (4.2.5) 
a) 


On the other hand, the deflection angle (8) is defined from equation 
(4.10) in reference [16]: 
(M_ -1) = 
8 = arctan | 2cot (8) a ren | (4.2.6) 
Mo {y,+cos (28) }+2 
The Mach number behind the conical shock wave (Mic) may, therefore, be 


- obtained using equation (4.7) in reference [16]: 
1 1 
a a I, 
Let M : 
: 1 2 n \ 
os ae i | eee 


nv M is a 
an oe 


The relation between the deflection angle (6) and the wave angle (8) for 
various Mach numbers are shown in figure 4.3 (reproduced from figure 4.2 


in reference [16]). 


4) 


guey _2¥ «gue _ 2 < g uez <9 ue} 
8509 (Ap = 1) = Oy  (oue4 T)o¥ = ,a¥ £6 ue} a a J.4 - a ale 
:Pauy a4nqdey YLUT 02 aALzeLay eauy MOLY 
AP[nuuy JaLUT JO UOLZeLND[e7) Aoy AyQauioay 2 "py aunbBry 


[91] suotynyos yo0yus anbiiqo e*py aanbL4 


[o/] STONY 3AM 
{09 s0€ 


wo 


.e) 


° 
O 
Do) 


(9) SIONV NOILO031330 


GV 


44 


A4.2.3 Area Ratio 
Using figure 4.2, the area ratio behind the conical shock 


wave can be obtained: 


A 
Hee) ey 


Equation (4.2.8) gives the area normal to the flow at the inlet lip. 
A4.3 Boundary Layer Loss 
A4.3.1 Mach Number 
The area ratio is related to the appropriate Mach number by 


by the formula: 


* Cy, )/2( -1)] 
A ryt2 | 
Ea (4.3.1) 
Le y.-l 
IC lis ae as 
2 1c 
Where A, is the area at the throat of the inlet. Similarly, 
(y.+1)/[2(y,-1 
At (,41)/2 ] Gt ifletys-1)] 
a Ee ee aa (4.3.2) 
1 Tm Seer a 
2.0 «1 


Where A 1 and My 4 relates to the area and to the Mach number ahead 


] 
of the normal shock wave. Dividing. these two equations dives: 


Be) 
dl 
o 


—i 
oO 


ae yale. 5 | rat C2041) 
= eure | (4,353) 
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Knowing Mico Ya: Aye Ay,» equation (4.3.3) can be used to calculate My i> 
indirectly, by subroutine CALM, which was mentioned previously. May 
should be supersonic (M4 > 1) to prevent unstart conditions. 

A4.3.2 Pressure 

The total pressure in front of the normal shock wave (Pr41) is 


received from the connection: 


Orin ay ale nD: eee 
Where Th» the boundary layer loss is assumed to be 0.93. 
A4.4 Normal Shock Loss 


The Mach number, behind the normal shock wave is defined as: 


2 ¥ 
jyye+et 
M aye — 
Me rae ee Serene t |, (4.4.1) 
12 | : | 
oi ale 1 


The total pressure behind the normal shock wave is defined as: 


Pra > Pra Seo 3 ye Coa - ie 
fart 24 / mal als! ea | (4 eees 
ea) oem 


A4.5 Subsonic Diffuser Recovery 
Similar to equation (4.3.3) one can obtain: 
| Heri os i 
ee Mine eae 
A fie \ate (4325) 
12 2 aa MR be 2 | 
2 v2 


Knowing May? Vas Ajos and Ao» the Mach number at the exit of the inlet 
can be computed using subroutine CALCM. The total pressure at this station 


is defined as: 


Pro = Pr12 Th (4.5.2) 
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Where the subsonic diffuser recovery (T)") is assumed to be 0.93. 
A4.6 Expansion Loss 

A4.6.1 Mach Number 

On their way to the combustor, the gases coming from the inlet 
expand at station 3; a sudden change in area from Ao to Az occurs. The 
sudden change in area acts as a flameholder by creating a hot recirculation 
region. In this section, the loss in total pressure due to this exnoansion 
is calculated. 


From the continuity equation: 


D_UoA, = pgUzA4 (1.4.2a) 
From perfect gas equation of state: 
eel! (1.4.3) 
PRT 


and from the definition of Mach number and speed of sound: 
M = U/a; a = vyRT (ied 5) 

Equation (1.4.2a) turns, therefore, into the form: 

P3 


2 RT. 


For sudden expansion of an incompressible fluid, the change in static 


Py _— ee 
RG My vYgRTo A 3 “Yakl3 Ag (4.6.1) 


pressure going from small area Ao to large area A. is given by (16, 25]: 


Ot Fon 


q 93 (1-A 


53) (4.6.2) 


Where Ao3 = Ag/ Ag. Also, for the incompressible case, the stagnation 
pressure ratio is given by [16, 25]: 


Pr3 - 5/Po 


ee ee 
iT? i+d9/P2 


(4.6.3) 


According to equation (4.6.3) as Mo decreases Pr3 approaches Pro: The 


model for sudden expansion with compressible flow is much more complicated. 
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It was assumed that static pressure is constant in the expansion. This 
assumption, which is reasonable for low values of Mo» is a conservative 


one, 1565 p=4/0-, 1S lower. 
Srai2 ‘/ 21 / os -| . D 
Substituting: T,/T3 = i + cary S74 5 2 


and p5 = pz, turns equation (46.41) G0: 


amo oe 
a 2 
po nS a eee (4.6.4) 
el 6; 
3 as ia M 2 
2 (lh 
Solving equation (4.6.4) for Ma: 
ae iJ 1+408-1 
oO. tae Ny 20 (4.6.5) 
Where: 2 
Yar! AoMs 2 


The subscript I shows that the computation was made from the inlet 
direction. 
A4.6.2 Pressure 


By assuming again that Po = D3; 


i 2 ie ee (4.6.7) 
Pro Pr2/Po _ 
and therefore: 
eae Deval ee) 
LAY aoe’ Soe eee 
Pr3r = Pro paren bee (4.6.8) 
web 


Where subscript I is as defined previously. 
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A4.7 Location of Normal Shock Wave 

The main problem in computing the Mach numbers and the total pres- 
sures at the inlet arises from the fact that the exact location of the 
normal shock wave is not known, and must be found. The way of solving 
this problem is as follows: 

First, solve for two extreme conditions by assuming that the normal 
Shock wave is located at the throat and at the exit of the inlet, respec- 
tively. After knowing the lower and the upper values for Mars Pr3r> 
iteration can be made to find the exact location of the normal shock 
wave. The criteria for this iteration is matching of values for Mo, 

Pr3 from both the inlet and the nozzle directions, i.e.: 


: _ (4.7.1) 
May = May» Pr37 = PT3n 
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APPENDIX B: TRAJECTORY EQUATIONS 

Bl. Atmospheric Functions 

Best fit curves were calculated for basic atmospheric functions, 
pressure, density, temperature and viscosity of air as a function of 
the flight altitude. The basic formula which was used for this process 
is as follows: 

F = A exp(-B X 107° n°) (1.1) 

Where h is the altitude in meters, and A,B,C are numerical parameters. 


The appropriate atmospheric functions are as follows: 


Py = 1.03322 X 10" exp (-59.148 x 107° hl-99) (7.2) 
oy = 1.228845 exp(-29.0144 x 1076 h |-15) (1.3) 
Tq = 288.16 expl-132232:K 1090 AE ae) 

When: h = 0 - 11,000m (1.4) 
Ty = 217.24°, when h = 11,000-32,000m Gum 
lg = 1-793 X 10°° exp(-45.1374 x 107° n?-°954) 

When: h = 0-11,000m (1.5) 

ug = 1.41724 X 107°, when h > 11,000 m (1.5a) 


In these formulae, the pressure (Dy) has dimensions of kg/m, the 
density (09) is given in kg/m, the temperature (To) is given in Ov. 


and the viscosity (up) is given in kg/(m-sec) (or: N«sec/m*). 


B2. Drag 
B2.1 Cowl Drag Coefficient 


2.1.1 It was found that the cowl drag coefficient has a strong 
influence on the results. Therefore, a new model for this cowl drag 
coefficient was developed. The model was based on a theoretical 


development done previously by Prof. T. H. Gawain [9]. The main 


50 


difference between this development and the classical theory, is that 
the boundary conditions are applied at the body surface rather than 
along the axis. 

2.1.2 The model, which originally was developed for simple cases 
(cones, etc.), was modified to fit the shape of the projectile, illustrated 
am Figure 1.1. 

2.1.3 The modified program is listed in Appendix G. In the combined 
Program (TRAJET) an interpulation procedure was used as a subroutine in 
order to simplify the calculation process. 

B2.2 Base Drag 

After checking the influence of the nozzle exit area (Ag) on the 
performance, it was decided to allow Ae to reach the maximum value 
possible (A.)> in order to reduce base drag. In this case, the base drag 
is negligible. 

B2.3 Skin Drag Coefficient 


The Reynolds number is well known to be: 


PgUyt 
= Ae AE 
Re, ; ( ) 
The transition Reynolds number of 
Re* = 2410° 


is uSually taken as criterion for transition between laminar flow 
(lower values) and turbulent flow (higher values). The incompressible 
laminar skin friction coefficient is related to the Reynolds number 


as follows: 


¢ = 1.328/vRe, 


DS,L (2 632) 


Where the subscripts DS,L stand for skin drag coefficient, and laminar 
flow, respectively. On the other hand, the turbulent skin friction 
coefficient (Cro 1) is calculated indirectly from the formula: 


op 


The computation is done in subroutine CALDC, which works in a similar 
way to subroutine CALCM which has been described earlier concerning the 
calculation of various Mach numbers. 
B2.4 Wing and Fin Drag Coefficients 
The wings and the fins of the projectile, also contribute to 
drag. Basically, each of these drag coefficients contains two parts: 
-Wing/fin wave drag. 
-Wing/fin friction drag. 
2.4.1 A psuedo 3-dimensional model was chosen to simulate the wave 


coefficient. The basic formulae used in this calculation are: 


ea MERE OUe one. Wee 
= et ee 7 S11) es ri = Pape Tae 
v(Mp) aes tan | ey (Mea-1), > otan J Mp zl ( ) 
a a 
y,*1 eV as ysl) 
Po/ Pro Se lat fe My” 2 (2.442) 
The effective span could be taken as (see fig. B2.1): 
fa & 
bi=b-3 (23453) 
Substituting: tanu = 2/c , Sinu = 1/M gives: 
b! Sy eS oC. teks 
y) Vier Dive (2.4 ay 
My -] 
The drag coefficient would, therefore, be: 
S052 (Pore = Poe! Bra Werb: 
COWL eee: Da Den D025 alae (2.4.5) 
iMag yoo! 10 Po’ PTO r 


The way this simulation uses the above ecuations could easily be 
understood when looking at the formulae together with the flow chart of 


the appropriate subroutine (Appendix C). Interference drag is ignored. 
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Figure B2.1 Schematic View of an Wing/Fin: 
b=span, c=chord, t=thickness 


2.4.2 The friction drag coefficient was calculated using the 
existing model for skin drag coefficient. (See Section B2.3) 


B2.5 Calculation of Drag 


Define the dynamic pressure (9) as follows: 


q= 4 ogy” (2.5.1) 


When q is in units of kg/(m.sec*), (or: N/m?) Also, define the following 
geometrical units: 


A = mR‘; S, = 27RL (25512) 


Where R,L are the radius and the length of the projectile, respectively. 
Similiarly, for the wings or the fins: 


Siw = nbc (2:55.53) 


Where n is the total number of wings/fins (a value of 8 was taken for n), 
and b,c are the span and the chord of the wing/fin (see fig. B2.1) 
Consequently, the drag (D) is given by: 


Dos ig J A2tC.. att Se. oC (C 


p “pn t Sp “ps * 5 


C (255;4) 


ww (Coww * Cowe? 


B2.6 Drag Coefficient of a Conventional Projectile Without 
Propulsion 


The program has an option to calculate also the trajectory of a 
projectile without a propulsion. The projectile is a conventional round. 
The formulae used to claculate the drag coefficients in this case are: 


2.6.1 Nose Drag 


Coy = (0.083 + 0. 096/M,*) (a/10)+- 89 
Where a is the cone half angle. 


2.6.2 Base Drag 


Cop = (0.6837 - 0.3165 M + 0.0525 Mo) (2/m) 
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2.6.3 Skin Drag 


Skin drag is calculated as discussed in section 82.3. 


2.6.4 Drag 
0D = Drag = q Ay (Con + Cop) + Sy Cos 


B3. Booster 

The projectile has an initial muzzle velocity of 2500 ft/sec. Part 
of the combustor volume can be used as a booster to accelerate the 
projectile even more so that starting the ramjet will be easier. 
Define exhaust velocity (U,) as: 


- 35.1 
Ve a: 9 


Where I is the specific impulse of the booster's fuel (in sec) and g 


sp,B 
is the acceleration of gravity (in m/sec). 
From Newton's law [6, p. 323] 


2 = dU 
F = th, U, * (m, - tit) oF (332) 


Where mM, and rhe are the mass of the projectile and the mass flow of 


the booster respectively. 


dt 
dU =nh, U -—-— (3.3) 
B “e mM, - hat 
-Consequently: 
m.- Mat 
AU = U(t) - U(O) = -U,2n —— (3.4) 
p 
Where t is the booster burn time. Hence, 
m 
AU = U, i, > Ty 0 (3.5) 


AU is the change in initial velocity due to the booster, where Mp is the 


mass of the booster. 


Oo 
om 


B4. Dynamics 


The flat earth trajectory with drag and thrust is well known. 


The differential equations of motion are: 


dnvaees ae 

aie q + (F-D) sin6/m, (4.1) 
m dx = (F-D) ) 

Pate | cos (4.2) 


where y is the altitude, t is the time, g is the acceleration of 
gravity, F is the thrust, D is the drag, 98 is the elevation angle, 
mM, is the projectile mass. 

For numerical solution of equations (4.1) and (4.2) a finite difference 


form can be used as follows: 


= - 2 - 
(F-D) cos At /M, - 2x54) X; (4.3) 


: 2 
[-g + (F-D) sine/m, ] At eye ay oY 4 (4.4) 


The values for j = 1 are from initial conditions, i.e. 


Kar anOs aVaie sc (4.5) 
Xo = Uy * cOS® * At, yo = Up * Sind * At (4.6) 
Veagoao Je 
g = arctan) #2 ~ Zi+1) = (4.7) 
j+2 7 *541 


Having calculated Xs40 and Vee from (eq. 4.3, and 4.4), one obtained new 
values for trajectory parameters using (eq. 4.7). Also thrust, drag 


and projectile mass are updated for j+2. 


Appendix C: Flow Chart of the Computer Program (TRAJET) 


€ 1. Main Program 
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choose: 


Right/ 
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NOZZ 
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120) 
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es 


choose: 
new Agr A 
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C2. Command Subroutines 
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| DIFFUS 


(Mo, 
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THROAT (2) 
(Mo. 1Pmo1 


EXPAN 
(M3r1Pngy 


choose: 


| THROAT (3) 
(Moy +Pmg1) 


| sR (3) | 


Piso” 982! 


| 


| DIFFUS (3) | 


EXPAN 


(M3q7Pn3r) 


Return | | 
| 
52h 
\ J 
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18% time ‘ 
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3. 


Individual Subroutines 
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| compute: 
| Png 


compute: 
gt 
TO 


67 


compute: 
hy sesieer 
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compute: 
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Thrust=0 


Return 
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(3.3. INLET 
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compute: 
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compute: 
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APPENDIX D: 
PROGRAM TRAJET: LISTING 
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APPENDIX E: COMPUTER PROGRAM LIST OF SYMBOLS 


PROGRAM EQUATIONS UNITS MEANING 
SYMBOL SYMBOL 


GEOMETRICAL SYMBOLS 
Z 


AREF A m Reference area 

Ad A, m@ Area at station j 

AIC Ale m@ Area behind a conical shock wave 
AJ1* As m@ Area ahead of a normal shock wave 
AJ2* Aso m ¢ Area behind a normal shock wave 
A30 Azo m ¢ Initial area at station 3 

ATAJ A,/A. - Area ratios 

L3 L. m Length of combustion chamber 
LPR Ly m Length of projectile 

RPR R, m Radius of projectile 

ALFA a deg. Inlet cone half angle 


*When: J=1,2.S, the shock wave is at station 1 2 real case, respectively. 
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PROGRAM EQUATIONS UNITS MEANING 


SYMBOL SYMBOL 
ATMOSPHERIC SYMBOLS 
TO Tg K Static temperature at altitude Y3 
PQ Po kg/m? Static pressure at altitude Y3 
RHOA 29 kg/m? Air density at altitude ¥3 
MUA Ug N.sec/m¢ Air viscosity at altitude Y3 
U0 Uy m/sec Projectile muzzle velocity 
MO My - Projectile initial mach number 
GRAV g m/sece Gravity (9.807) 
GA - - Air heat capacities ratio (c,/cy) 
GA - = (y,+1)/2 
GA2 - : (y,-1)/2 
GA12 : (yt )/[2(y,-1)] 
GA3 - - ee 
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PROGRAM 
SYMBOL 


RHOF 
ETAT 


GF 
GF2 
GF12 
GF3 


EQUATIONS 
SYMBOL 


UNITS 


kg/m° 


MEANING 


COMBUSTION CHAMBER'S SYMBOLS 


Fuel density 

Burning efficiency 
Burning rate parameters 
Air mass flow 

Fuel mass flow 

Total mass flow 

We/W. 


Burning rate 


¥34Re*T Ty / T(when:T = T(y¢)). 
: fe) 
R(j/mole/ K) 


Hot gas constant 2 
MW(kg/mole)*g(m/sec 
Hot gas heat capacities ratio (c,/c,) 


(y_4] )/2 
(y_-1)/2 
(ye+1 )/[2(y~-1)] 


Yel (Y¥e-1) 
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PROGRAM 
SYMBOL 


TJ 
TTJ 
PJ 
P TJ 
MJ 


PTJ1, PTJ2 


MJ1 ,MJ2 


T3M 
M3N, M3I 


EQUATIONS 


SYMBOL 


UNITS 


kg/m 


MEANING 


THERMODYNAMIC SYMBOLS 


Static temperature at station j 
Total temperature at station j 
Static pressure at station j 
Total pressure at station j 


Mach number at station j 


As above with AJl, AJ2 


Maximum 13 available 


Mo calculated from nozzle and inlet 


direction, respectively. 
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PROGRAM 
SYMBOL 


EQUATIONS 


SYMBOL 


UNITS MEANING 


LOSSES SYMBOLS 
- Boundary layer losses 
- Subsonic diffuser recovery 
- Nozzle losses 
- Conical wave losses 
- Normal shock losses 
- Expansion losses 


- Heat losses 
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PROGRAM 
SYMBOL 


CF 


Thrust 


ISP 


SFC 


EQUATIONS 


SYMBOL 


orc 


UNITS MEANING 


RAMJET PERFORMANCE SYMBOLS 


~ Thrust coefficient 
N Thrust 
(or kg) 


N/kg.sec Fuel specific impulse 
(or sec) 


kg/hour/N Specific fuel consumption 


Lis 


PROGRAM 
SYMBOL 


Drag 
WPR 
WB 
DELU 


EQUATIONS 


SYMBOL 


UNITS 


MEANING 


TRAJECTORY SYMBOLS 


N 
kg 
kg 


m/sec 


m/sec 
sec 
sec 
sec 
deg. 
deg. 


Drag 
Mass of projectile 
Mass of booster 


Change in initial velocity due to 
booster 


Uy + AU 

Change in time 

Time 

Time of burning 

Gun elevation angle 


Projectile second cowl angle 


DRAG COEFFICIENT 


Nose drag coefficient 

Wing wave drag coefficient 

Wing friction drag coefficient 

Skin drag coefficient (laminar/turbulent) 


Base drag coefficient 
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PROGRAM EQUATIONS UNITS MEANING 
SYMBOL SYMBOL 


MATHEMATICAL SYMBOLS 
PI TT ~ 3.14159 
IPR - - Printing parameter: > 0 combustion 


results together with trajectory 
(on different files): 


ul 


0 working results only; 
= ] also reasons for not running 
= 2 also full reasons for not running 
= 3 also 100p on mach number (CALCM) 
= -] trajectory prints only 

ITRA - _ 2 Loop paramerer: 


= +] single value for Ao/An» Ac/A., 


-] loop on Ao/ Apne Ac/A.» and 


print summary, only. 
ILOO - - Check parameter: 

< ] regular run 

> ] doesn't have a solution 
fin 11.0 - - Trajectory printing parameter 

(prints every ILO point). 
TRAM - - Ramjet parameter: 


= 0 ramjet in operation 


1 projectile without propulsion 


XMO Kuo - Stopping mach number 


bs 


PROGRAM 
SYMBOL 


INIT 
BURN 
TERMO 
INTER 


NOZZ 
CALCM 
RESUL 
CHOKE 
INLET 
CONE 
THROAT 
NSR 
DIFFUS 
MIN 


~CORVAL 


EXPAN 


HEAT 


CHECK 


TRAJ 
ATM 


EQUATIONS 


SYMBOL 


MEANING 


SUBROUTINES 
Computes initial conditions 
Computes combustion chamber's performance 
Thermodynamic comand subroutine 
Computes, by interpulation, thermodynaic 
conditions, or cowl drag coefficients. 
Computes nozzle performance 
Computes mach number. indirectly 
Computes ramjet performance 
Checks if nozzle is choked 
Command subroutine for inlet 
Computes conical wave loss 
Computes boundary layer loss 
Computes normal shock loss 
Computes subsonic diffuser performance 
Command subroutine to compute situation 
when normal shock is at point 2 
Command subroutine to compute situation when 
normal shock is at the correct place 
Computes losses due to expansion into the 
combustion chamber 
Computes heat losses at combustion chamber 
Check for pressure capability 
Command subroutine to compute trajectory 
Computes atmospheric conditions as a function 


of altitude (y) 
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PROGRAM 
SYMBOL 


BOOS 
DRAGG 
CALCD 
DRAGN 
DRAGWW 


PRANT 


APRANT 


PRES 
DRAGPR 


DYNA 
DATTA 
PRILOO 
PRIN 


EQUATIONS 


SYMBOL 


MEANING 


Computes booster performance 

Computes drag 

Computes skin drag coefficient, indirectly 
Cowl drag command subroutine 

Computes wing/fin wave drag coefficient 
(command subroutine) 

Computes Prandtl-Meyer angle (v) 

from a given mach number 

Computes mach number from a given Prandtl- 
Meyer angle (v) 

Pressure ratios formula 

Computes drag coefficients of a projectile 
without combustion 

Computes the dynamics of the projectile 
Initial data 

Prints ramjet performance 

Prints detailed values when does not fine 


solution 
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APPENDIX F: COMPUTER PROGRAM USERS GUIDE 

Input Data 
AOAR, A5AR, A1AOQ, A2A0, A3AR 
ALFA, TEA, Tel 
IPR, ITRA, XMO, IRAM 
PID], PID2, PIN 
Options 
ALFA = Inlet cone half angle 
TETA = Gun elevation angle 
TETP = Second cowl angle 
IRP = 0 clean print of RAM + TRAJ 

= 1, 2, 3 more details on RAM 

= -] TRAJ only 
ITRA = 1 works on one set of data 

= -] loop on AO2A5 
IRAM = 0 Ramjet 

= ] Projectile without propulsion 
Execution Commands (For use with IBM 370) 


F2.1 Opening Commands 


LAXXXXP 
Password 


GLOBAL TXTLIB FORTMOD2 MOD2EEH 

F2.2 Compilation 

FORTGI TRAJET 

F2.3 Run on Terminal 

FILEDEF O02 DISK TRJ D(RECFM F BLOCK 132 PERM 
FILEDEF 03 DISK DRG OD(RECFM F BLOCK 132 PERM 
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PILEDEF O05 DISK INP 0D 

FILEDEF 06 DISK CMB D 

EXEC RUN TRAJET 

XEDIT CMB D 

: XEDIT TRJ 0D) 

: XEDIT DRG D) 

Note: CMB D will contain the combustion process results. 
TRJ_D will contain the trajectory part. 
ORG_D will contain drag coefficients. 

F2.4 Hard Copy 

FILE 

PRINT CMB D 

PRINT TRJ D 

PRINT DRG OD 
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APPENDIX G: 
Gl. PROGRAM acRo'*) : LISTING 


(*) The original program was developed by T. A. Gawain [9]. The modification 


listed here was prepared for this report to calculate the cowl drag coefficient. 
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ia) aA aAANg 


eo a) 


ILE: AERN FORTRAN A NAVAL POSTGRADUATE SCHIIL 
AF200010 
SPECIFICATION STATFMONTS AF29 020 
te *~ 4£99 9930 
IMPLICIT RFAL®4( A-H,0-Z), INTEGERS 4(I =") AF2N90460 
API9NNQ059 
OUMENSTON X( 702) 491792) 4RP(7O2 XE 702)9FO1702) .GP( 702), 4F290060 
#6(702) 5 QA0(702 340420792) ,ONALTO2) pOMA( 707) KAD (B) 4£99°1970 
AFT SCHOO 
DATA KOE SONNE t/ KT Le AEA tafe ee SA UA Al Ba AFP 99090 
FKPAOSEPIESKOTS FAT US SKINS TONES, KNT/EOQ TOS gKBPA CIVIL MUST SOL AFO79100 
+ eT ANAT 2037 2NK/ 200) » NTAB/ 207 sKNOS EFEONE, KTATESIOTUS, ACQN4aLLo 
+#KAP/*X CC® ,PARCTY OMATES, "(TOS I NEGSH, TTI VSt,t) ‘7 AEVOOL2D. 
+e DL /=0 00 7300/ 0 CL 2/0 cl 7707 NCJ bv, NMB2/OL/,INASES OF AFTINNLIO 
AFONIL4O 
SELECT OFTIONS ATQ)W1SO 
AG?9 160 
1001 WRITE(6 21020) AE2VQL7O 
1021 FORMAT(C PVE," ENT=HO CANTRAL CODE AS FOLLWS:*7 AF29) L949 
+T6,!'7T9 set hin EX°' ANAT ION  T4L,tENTER HENS AFI2N %199 
+T6e9TT ENTER FNPUT PST AMETEPSI,TGL,ZIFNTER INOS ED() Y299 
+TOe' TS CLO4PUTE GUTMETRICAL ARQAYS',T4G1,*ENTEQ CGI/ ACID I20) 
#T6e9TA CALCULATE TUTDIT® "OL," FENTTR CIS AZ=Q99 7220 
#T6o'T7 TAHIAATE OUTPUT, ToL, INTER The / AERN 729 
#TGe* TA PLOT MUTONIT® ,T4aLet "NTE OVE AE™7 0240 
#T6e'TT QUIT PROGRAM! ,T4L,"ENTER QTY/ ) Az2709250 
PEAN Se» LC4L,FRPQ=1NAL) KAR AFO9)200 
1041 FROMEAT(A2) AFO07979 
TFE(KAR SE CKIF) G TN 2001 4529780 
TF (KAP SFC.KIMD GP oTN 208) AE7% 9299 
TF (K29 fC 2405) GS TO 274t 4£999300 
TF (KA9 ,FC.KCO) GE TN 34281 AE®9 3340 
IF(K@R FC .KTN) GE TA 33A1 AF?) 9320 
[F(K4R.FC.KPN) Go TN 3601 . AEI')339 
TE(KAR CEQ BKKIOT) Ge TO 3721 AEQ) 1349 
LO6L WRITFI6,1081 ) AF2.99350 
LOBEL FORMAT(® #,) INPUT FR2N0R, CEPEAT INPUT?) AEPQIHY 
GOT? Lucl Ree one 
AF29)9380 
EXPLANATION 4£990923290 
AF70 0400 
200L WRITE (6,2021) AZ2NNG10 
2021 FURMAT( ® #'.* SUPERSONIC FLOW ABNUT A SLIGHTLY Tae oe SLENNER, '/ AFD!) 1420 
et @ PONTATEN BCNY NE CEA ahead THIS PROGRAM UTILIZES Ats A£9\) 5420 
+8 Pt, OMNNTEIS LIM -ARIZFN ALYSTIS DEVELOPED KY T He GAWAIN, 8/ NFO) 440 
#* *,9ANA MONTFIEr FOR ChwWL NARAG CALCULATION DF A GUTONE DY TO 9459 
ee UC APQAOJECTILE EV OcAMICHAT .* 7 AEQN 4460 
8 "SCL ASSICAL LING ARITZED THONRY SATISFIES Hoe BIINA AOE? AEPN 1470 
© H,ICONCITIOCNS CNLY APPROXTMATFL Y, AT THE AXIS: THF SIFSENTIY AF )94A0 
+8 §,9ME THON SATISFIES THEM VERY ACCURATELY AT THE APDY SYRFACE'Z AFIONN49N 
#! 8, °*rTS€uFe.*/ AFQ9 9500 
+?) TRE BONY IS SCALFD SO THAT THF LAFGFEST CROSS=$SFCTIIN IS#/ AE999519 
+8 "598A CIRCLE OF UNTT RADTUS. AERDOYNAMIC COFFFICIENTS ARE BASENP/AFIDNG ID 
o* 8,9 FRONTAL \QFA,TS/ AFD 9520 
+! oye TEE CONFIGURATION CONSISTS OF A CONTF AL OR AWIVAL NOSE, */ 47299154) 
#6 F,080FNCTTN SY “Y*¥AOLS CN OR ON, RESPECTIVELY, SNES AEST 550 
et! USINOTINNAL MIPRSECTION TF UNTFOTM INET SAOTHS , ANA AM QOTIONMRF/AFS 111560 
©? @©,0Af TAIL IN TRE FORY AF A TRYNCATAREN CONE 18 TWEE, Sy4MANT FTIVALI9D570 
#2 6,°9NA CT, GFESPECTIVELY, TOTAL LENGTH XL TS nya pwraen TNTO NL OZ AFRAI5AN 
#8 OC SHENTEP VALS OF OUINTFORY WINTH 4=KXL JL. NOSE LENGTH ANON wT ne]ts 4£9;)C890 
e+) OPS CC TING PT ENGTH ADE EXPRESSED AS MAEM AND NABH, I COECTHI/ AF2% 600 
#* 8,8 VELY, WHERE NA AND “!H_ ARF INTEGERS. PARAMETER CL DENOTES!) AE~0 9610 
RITS(6,2041) AE2NN620 
2041 FORMAT(® *,* THE (eCRSASE IN RADIUS JVEP THE TAIL Sree sf AERN 620 
+t 9,! CALCIAT [MNS ARF “ADE FO EVETY STATION J BUT OUTIUTYY AFI 3640 
+" 8.9TS TARULATEN CMLY £N® Ja2eK®NTAR WHERE KzVe loo devee!/ AF20 1650 
+f ote MACH NUMP-P 4 MUST LIF WITHIN THE OANGE Te ceMcty AE RN 1660 
#8 FLPSORTI LS. 2 OMAX€#%=2) WHERE © PMAX DENOTES THE MAXEMI AINYES AF70 1076 
+! ty tery ope,ty AE 299640 
ities! INPUT DARAME TERS ARE XMy XL, Ole NLe NA. NBy MTABLCN ORP/ AVEO I4N9 
+? ©,10OM, CT OR OF, AND AM OPTICMAL ALPHAMIMERTAC [OEMTIA=9/ ARDY YT 
+8 (SPE FER 1D TA a CHARACTERS TN LENGTHY / ASQ) GITY 
oo aye A CONTENM VOINT IS LOCATE) AT KXe0ee AT THE ALNPOLENT OF TZ AFIQ 20 
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SRADUATE SCHOIL 


NAVAL POST 


A 


FORTRAN 


AERO. 
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G3. PLOT ROUTINES FOR USE WITH AERO [9] 
-PREPLOT P (PLOT ON PRINTER) 
-PREPLOT G (PLOT ON PLOTTER) 
-CHARTS (CONTROL) 


Loe 


FE 


a a Anann™s 


TLEs 


7AM 
7 EXE 
/ FOR 


~ 
° 
°o 
o 


teers JINB ¢ 


PREPLOTP FORTRAN A NAVAL POSTGRADUATE SCHOOL 


2250+0258),'AMICHAI X2936' »CLASS=A 
FRTXCLGP 
Te Csysis co * 
SPECIFICATION STATEMENTS 
IMPLICIT REAL ®4(A—H»0-X) ,INTEGER*4(I=N) 
OIMENSTON X€502),R(502)5 QNA( 502) 5 9A0(502) »QA2( 502) »KAP(8) 


READN(5+1CO0) XMeXtLeOLeNL »NAo AB pNTAS, KNOSESKTAIL, 
+KAP,CNA,ZCAQ,CA2,X ac 

FORMAT( 1% +35 14.6 5640 2X0 13) 620 2K AZ ISL Xe BAG/LXy 4E14.5) 

NLP2=NL +2 

50 1040 J=lyeNLP2 

REAN( 5.10290) KIS) RID) QNAC I) 4 QA0 (J), QA2(J) 
FCOMAT(LX -SE1L4.5) 

CONTINUE 

wQITF(6,1060) 

FORMAT(I'1°,* SUMMARY OF RESULTS*///) 

WITTE(G 41080) XMeXL yO sNLoNAwNBeNTAB «KNOSE,KTAIL » KAP (6) ,KAP(7) »© 
#CNA 6 404C42, KAC 

FORMATO PC XMS FT Se TILA ey  XLa fof eG TIO oe WN! oF 60% Rage eee 
#'O',*NaAa’,[3,TL8e'NB=!? 1347369 'NTABD=!',13+T5S4e'NOSE topes: ’ 
#70", "TALL CONEH",A2,T3b, 1D=", 2A4/°O", 'CNA=" oF 6s 4» Tie escages, 
FFE GeT 3G 6 CAL= "oF 7a ool 54 e' XALE FF HL3) 

aRITE(6,1100) XM,KAP(6), KAP 7) 

FORMATUEL 36x ¢'SADIUS & VS XK XME* GET SGy! 105", 2A4///) 
CALL SLOTP(K »28 »NLP2, 4) 

WRITF(6,1120) 

FORMAT('1',3 0X ,"*NCRMAL FORCE COEFFICIENT QNA VS X*/) 
ARITF(641140) CNA,XAC,KAP(6) sKAP(T) 

FORMAT('C*,36OXetCNA=", Flo ay? XAC=*,FI.G,?® 102" ,2A4¢//) 
CALL PLOTP(X ,QNAeNLP264) 

aQ@ITE( 6.1160) 

FORMAT(C P12", 36Xe"AXTAL FORCE COEFFICIENT QAO VS X*/) 
WRITE(O¢1190) CAQ,KAP(6) sKAP(T) 

FURMAT(9C',30X9'CAQ=', F743! 102 ,2A4//) 

CALL PLOTPIX ,QAD sNLP2, 4) 

wRITE( 6.1209) 

FORMAT (914 ,36K,*AXTAL FORCE Coe rae ENT QA2 VS X'/) 

WRITE(5 41229) CA2.K A916) KAP (T ) X 
FORMAT (*C*,36X,'CA2=",FT eae! 1O=* ,2A4//) 

CALL PLOTPI(X,QA2,NLP75 4) 

STOP : 

END 


/* 
//4GOeSYSIN DO * 


133 


vv 
pv 
mn 

folelelerele) 


vuyv 
DRA 
COOGOSCOSCSOCHOSCSOSSOOSCCOSGGOSCGSICOCOOOCCOSOCoSCOCOSSCOeC0Oo 


Vv 
mR 
TANT NAN AM ADNAN AATAARADAMAMADNH 


SCOGCOCOOSSSCOCOCVOSOCOGSOCSOSCOTGOOS CE OVGOOOOSCCOD0O 


vv 
2D 
COCCSOCOGSCOOCOCOGUOCEESCOGOCUOCOGSOOCCOCOOCOSOSG9 


vv 


PR 


PQ 
PR 
aQ 
PR 


evvnvue 
RDOVDODXDAI 


evuv 
RBRODBAD 


vvuvev”e 
DDODDRAD 
anannn 


v 
R 
m 

DANONE WN RH DOONAN EWN O DBWKNEWIVKODDIRNEW 


FILE: PREPLOTG FORTRAN A NAVAL PQSTGRADUATE SCHOOL 


4AM1C3250 JNB ieee or 0258),'AMICHAI X2936* ,CLASS=A PREOOOLO 
7 EXEC FRYIXCLG PREQ0020 
/FORT.SYSIN gS * : PREQ 


SPECIFICATION STATEMENTS PREO 
IMPLICTT RFAL®4(A-H,1—X) p INTEGER®G( I =N) PREQ 
DIMENSION X€ 502) 421502). QNA( 502), 940( 502) »QA2(502) +KAP(8) eRe 
RFAN(5.1CC0) XMeXt OL NL »NA,NB»NT AB, KNOSE,KTAIL, PRE 
+K AP sCNAy EAQs Gaz KAC . PRE 
1000 FORMAT (1 Xs 3E14 504 2Xo 13) 20 2k eAZI/1 Xe BA/ 1X 4El 405) PPF 


2 
J 
2 
T 


a a 


) 
0) 
9 
D) 
fe) 
0 
9 
Y X(S)eR(D de QNAC I) ¢ GAD (J), QA2(U) PREQ 
e15 5) PREO 
PREQ 

PREG 

SUMMAOPY OF RESULTS'///) ppcg 
My XL» IWeNbLeNAeNBoNTA SEeKTAIL »KAP(6) pKAP(T) 9 pees 
0 

fe) 

D) 

0) 

9 

4 

a) 

) 

0 

Qo 

9 

) 


-“~- th 


co eee 


1020 
1040 


1060 


~ “pO -& 


2ON0eee 


eOL=® pF6en%sT5Se'NL=* pi 3/ PRE: 

OAS Blab CODE=",A2/ PRE 

NA=",FGe%sT18,'CA0=", pRe 
PRE 

R 

3 

3 

3 


Be2Ze 


27 2 = PONPTN-D 
eo + 


VPUVUVUAA Zap 2 


‘ 

’ PRE 
9250050500. 98234 -) ORE 
’ 


DRE 
PRE 
i 
1760» SYSIN O00 #* PRE 


8 
4 
= 
0 
3 
A 
, 
‘ 
4 


ROO0~ < -- 


HOCH WP PmODeKanc £90 
DVUUW =“N == 


ODVS9D0— tl 22 « ee TD 
© XK KM DOW = rs) “Oo *O- 
Gee 22 NO KK~Y~ 


Omer Rene! NUN 
<« COOVUNAED 


x 
Cc 
= 
L 
' 
’ 
NX 
A 
4 
0 


ee ae 


VErOOO es + Fo 4SAAND 
KOON - 


134 


DOOHO FT OOLT FOTOS OF DODIINONDO FNIFSOAIODNTS RADA LPHOAAHR 
ANDO 9 472 DA HOS DNHA DTT ITI ID Ft tS SIDVO SONI 30=fE990 
C0900 9292979003999 95 99392999095039595350 70399599 
00090904 299 999C38 979999E9279999SS99C 9030f29990C992 
OIDANDNIMI™ DIRDANONTN OF VO DHNDNFNOCHDFSANMDFNOEDND 
WVO09 DCD DDD tt 4 A IIN ISIN NANI OUI II 
rYrOVVIO DV§I VII VOX SV OVV VV Ss9S CCV TO S.OCC so XONCSeSVS 
026909929290 9591290990997T292T99SSC9N9F 3990090090990 
++ > + 
osoocsa 
ANITA » aad 
o20n"7 r a J aq 
os339 —“ VY Alva) 
oo9000 + =e >> ~ 
Wu uu - oO ale) 9 
rxrrxre ensjo= "y a 
aaaaa So ea be Qa 
« wIdOr = =) 
La AW z= w 
vw uwye@ - - <a aeaazreuw 22> - 
~ or - - a Ls >i 7. > 
tO " vy ce ete — . 
- mn fon} “vn > ea) -— - 
< nne 2 aw NNMAYA I bog a | “ 
x wut a " > rm oe > 
- NNO OD - Fe Vale eaalt a) mw ww s 
a 1 DO ~ ~- viehin~ =e LJ 
wu uw> beg h Z& aeeZZ ® Se z 
ae “uw wv a2 nunw4A— “es 7) 
a oP ~~ @ “ e Ouod04N xm uo 
" co omer uw a sss ee feu rU = 
4 Wo Him W om is Dbalatar 4 ——_ ipa) 
“ve va on2 = >» 39 «x osoo4 wi i=] 
au .= 1S E> am oN Lr octcevoo.s ww o 
+ ex age Vv x4 ND NNMANE [Ss) a) an 
vu ar eon WwW ~ . " >P>>oHy ad > wo 
vo = a~a~= oo a FETE aa =z Londaed 
-« Of =--9O oN A Wn ww nn " (STS) 
ow NN? SO oe ee A: 4° I~) o> t 4 Wu 
<q am Oo eee One ST ~~ WO WWW, -—_ w >> 
a oO zw AND H et eee NMA aw a) wu 
“ N We ~——-y J «© 24 a Lana HHnngo wn Piaaed || IW 
“ww x FS ™M“ ee YY WIN WO Se Pe 4 aa Pe nu 
a xO S sayr wu OND Wo YUUOUY aq aU Zz 
a — ww vt >>~O Xt = Om Ot S>>>>FYU ee - un 
Vv ! <a < 2 VOUT jor -4~ 2 2.8 ae = Tre ou 
= x rw " ~~o “Dy -—ut D> Yerro wu ~a ae 
° YW & min One Dw ww LItLw au, vn -—_ 
o = 292 u Wu er up DY QOD YMUvwesa wn wQyv nw 
Oo x Zz a. VOVUH-ON~ FIND HH Hho “eH Enzi uve 
o < a o aAqtwue Oo 4 aaaan [20 efgwoa wmMuw 
NN e z w eouwIimg rum AU YWNYVIN>S ee ANNOY WY 
x - ar uc Mviaqaaquwte ~~ 8 NU LZR =—N Qe uw 
i -~ aon we = A aANUY eee DONO ow ze -M10 aa 
~ o -O ac AN HKD wt WOIA Dee & eM ee Ue & eee 
< wo aw wt Qo roMNn es SOD © > cCODCwW> VU 2aadao a00 
x N OF Wat NMAIIYL DMI FOMM LAIN NWN KNEE OWUFBTBEN Fwy 
Vv Oo ux aw KS KUM2UV CHF KLMTMINMMMMIUNNH HN > OHH Hem NOU 
~ e r 4 —ur WAI —AC wr HYIY) DL Pwr mm aie OU BOR WOU Ie ee 
=z Oo Wee SIH POH DUH IDYOOTH HHH eLIDVIDZ2ONKH DDH Hy 
a Ww LY UONnOzrrrHaNnsHOce KH PH taker HH CLCHHOAMONrOUOCLaADZ 
. N Fe OFNV HK AND FW VOSA FS F ZV ZZAMZS® Phe NNNE, 
i? MA © CAKYDZZVHY OFA DVNVINNDIZZZZOSCH HNN HOF OZ>Y Ore da 
~ ~OMNW AAWVODIOVNO0A MWADDVODOONDIIDIAUNYINIOMDO AIMMOOO 
o as [ax | 
w ou - a 
N OoOxun vV900000 VO0 A BOX0D0DOD0VO00R0O0O0000U0009000 
°o Ie lalelalale) ZWOD 9 ANDOOVOCWOHNDOOAD waqgqnaa 
e “ra ad kat < ~< 
° oOurFv Wwe WU be Weta azw ~ w 
ww « w Oo £2 ~ 20 o°00 oo~ 9AOrF a 
N No 4a SM LINZ YO =~ 2 oO COODANOCYH —OVUHNA 
Leals 9 aenovu WD Wha *. ect ft uw — WUUKYIr wor 
~ou UwU xrarsvs RFuus YY ~ WU rr RN Or 
pare — =O Kavvniu awuyuy ve ww “” OVOVYLY *MUWUYIUUY 
mu waw auU>r>>r VUtrr § FF Py Ur were Sr 
Ux wet ee euuviviwiviws uwvuevws vr ve ve Se) A ee ee 
Tre NN RF ON RO RFE OK OE OE FON RS 
ao~ NN OR OO ON ME OE CN OC OC 
ouv £ 
“wo >- = 
nuwn 
mu * > 
owe ; 
eo oo NTO DOO OANM FF NOFDKOKNMINOFDTOAKNGATH OF 
o 


4/ 
4 
4/ 
/*® 
4/ 


135 


A 


EXEC 


CHARTS 


FILE: 


NAVAL POSTGRADUATE SCHOOL 


FILE FTOLFOOL A PLOT FORTRAN A (REPLACE 


! e 
aie [4 
w 
| oo 
(‘= 
Oo 
a 
arn 
ww 
ero 
Ze-O 
Oo 
cur 
aaw 
tte 
Zz 
cc 
Fre az 
coo <a 
oo za 
uwu> ar 
uw (<a. 4 
oa -o 


Odserg wWhiwnooud 
aww be - om 
ee 2 ee SWOr wou 

Zz QAw-FrOT >U 
Muir xadouvdauaw 
War b WA ZNO ZOO K 
WII LI WIET | ILD | w | OW 


136 


G4. Results from AERO/COWL 

G4.1 The symbols used in programs AERO and COWL are defined in AERO 
and are presented in Figure G4.1. The values of a, bys Cy; Do» Cy are 
normalized with respect to rm. NA, NB, NC, NB.» NC, are the appropriate 
numbers of points used in the program. The cowl angles GQ) > %, were 
selected as 20°, 25. respectively. 

To create flow at the first cowl in Figure G4.1 which is the same as 
for a ramjet inlet, an extension to the body was used. The extension con- 
sists of the cone of length a and cylinder of length by in Figure G4.1. 
The cone angle, By> is 8 degrees. The value of b, was varied until the 
pressure at the first cowl was equal to ambient; a value of b, equal to 
vr) gives this condition. 

G4.2 The normalized values of the various variables which were 


selected are as follow: 
az/7.11 , by =e Cy = O.2: s bo = 2.104 Co = 1.06 units 
- DL = 0.0736, - OL. = 0.177 units. 
The appropriate numbers of points are: 


= 61, NC, = 30 


NA = 203, NB = 200, NC = 6, NB D 


z 
The dimensional values are: 


Cc, = 0.40, b, = 4.22, Co = 2.1 inch 


2 
2.48 inch 


5 
4 
iT} 


1.985 rp A r3 


G4.3 It was found that the cowl drag coefficient is sensitive to the 
magnitude of projected cowl area. For example, another combination of 


these variables: 


L3/ 


aS Wilt b, nt he Gai= 0.38, b,. = 2.16;%¢ 


9 = 1.46 


2 


NA = 196, NB 


193, NE .=ah1, NB. = 60, NC, = 40 


=DL = 0. 138(20°), =D. = 02283. (9.5-) 


2 


Gives higher values of the cowl drag coefficient: 


Mo 3.0 223 


0.0732 


0.0953 


0.1214 0.1562 


Therefore, attention was made to select an inlet shape as smooth as 
possible. 
G4.4 AERO can produce also graphical results. A typical example is 
illustrated in Figure G4.2. 
G4.5 Computer Program Users Guide 
Use opening commands and compilation similar to that defined in 
FQdget eee. 
The routine to run AERO on the terminal is defined in program 
itself. 


The routine to run COWL on the terminal is as follows: 


FILEDEF 02 DISK OUT D (RECFM F BL@CK 132 PERM 
EXEC RUN COWL 
PRINT OUT D 


t optional 
XEDIT OUT D 
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H1. SUMMARY 
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0.25-0.40 


0.25-0.40 


0.25-0.40 


0.25-0.40 


0.25-0.40 


0.25-0.40 


0.25-0.40 


0.25-0.40 


0.25-0.40 


0.47 


APPENDIX H: RESULTS 


Ao/Ag A/A ) 


r p 
0.887 0.426 oso 
0.827 0.426 9.9 
0.827 0.426 = 
0.75 0.426 935 
0.827 
0.91 
0.827 O27 gn 
0.32 
0.47 
0.827 0.426 6.5 
pe ee) 
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